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The optical analogues of Bloch oscillations and their associated Wannier-Stark 
ladders have been recently analyzed. In this paper we propose an elastic realization 
of these ladders, employing for this purpose the torsional vibrations of specially 
designed one-dimensional elastic systems. We have measured, for the first time, 
the ladder wave amplitudes, which are not directly accessible either in the quantum 
mechanical or optical cases. The wave amplitudes are spatially localized and coincide 
rather well with theoretically predicted amplitudes. The rods we analyze can be used 
to localize different frequencies in different parts of the elastic systems and viceversa. 

PACS numbers: 43.35.+d,63.20.Pw,43.40.Cw 

Recently, undulatory systems showing analogues of Bloch oscillations and Wannier-Stark 
ladders (WSL) attracted increasing attention in several fields of physics jll 0, [1 U, Ej]. As 
shown by Bloch, electrons in a periodic potential have extended solutions. The same is true 
for the behavior of an electron under the action of a static electric field. In contrast, and 
opposite to intuition, when both the periodic potential and the electric field are present, the 
solutions are localized; this is only true when band to band Zener tunneling is negligible or 
the system is short enough. The spectrum then shows equally spaced resonances known as 
Wannier-Stark ladders, the nearest-neighbor level spacing being proportional to the intensity 
of the external field [6|. In the time domain, the Wannier-Stark ladders yield the so called 
Bloch oscillations which consist in a counterintuitive effect where the electrons show an 



*Permanent address: Institute de Fisica, Universidad Nacional Autonoma de Mexico, P.O. Box 20-364, 
01000 Mexico, D. F., Mexico 



oscillatory movement under the action of the static external electric field 



ELQ. 



The predictions by Bloch and Wannier gave rise to a long controversy that lasted more 
than 60 years for the Bloch oscillations |9| , and more than 20 years for the WSL Q,HQ. 
The ladders were observed before Bloch oscillations. This was done first in numerical ex- 
periments using simple one-dimensional models k| and later in the laboratory Bloch 
oscillations were also observed later on |l5|. The most important ingredient to explain WSL 
is the wavelike behavior of the electrons. Therefore, these ladders could also be observed in 



classical undulatory systems. Some of these classical systems have been analyzed theoreti- 



cally 
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In this paper we study special elastic rods whose torsional waves for free ends present 
some analogies to the WSL. The first system, depicted in Fig. (a) and which will be 
referred to as system A, consists of a set of N circular cylinders of radius R and varying 
length l n , n = 1, 2, . . . , N, separated by very small cylinders of length e <C l n and radius 
r < R <C l n - This is the elastic analogue of the optical system with varying widths used in 
Ref. [l|. System B, shown in Fig. ^ (b), is a beam formed by N cuboids of constant width 
w and length I. They have different heights h n , for n — 1, 2, . . . , N, with w, h n <C I. These 
cuboids are separated by small cuboids of dimensions h' = w', e' < L This is the elastic 
analogue of the optical systems with a gradient of the refractive index along the direction of 
propagation used in Refs. 0,0,0], as we shall see below. Systems A and B were constructed 
by machining a solid aluminum piece. 

We first discuss the design of these systems and then, from a qualitative point of view, 
the normal-mode frequencies and wave amplitudes for torsional vibrations. We later use the 
transfer matrix method and obtain the normal-mode properties, which will be compared to 
experimental measurements. 

In order to design systems A and B we start with what could be called an independent 
rod model in which each body oscillates independently from the rest. The normal-mode 
torsional frequencies f ■ of rod n with length l n and wave velocity c n are given by the 
well-known expression (19( 

ff = (1) 
where j is the number of nodes in the wave amplitude. To get a set of equidistant frequencies 
we consider two options: either the lengths l n are varied with a fixed wave velocity or the 
lengths are kept constant and the wave velocity is changed. For system A we take circular 
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FIG. 1: Rods used to obtain the Wannier-Stark ladders: (a) rod with varying length cells and (b) 
beam with equal-length cells but different heights. For system A, l n = 1/(1 + 717), n = 1, .. . , 14, 



with I = 10.8 cm, 7 = 0.091, e = 2.52 mm and y/G/p = 3104.7 m/s. The radii of the small 
and big cylinders are r = 2.415 mm and R = 6.425 mm, respectively. In system B, I = 5.0 cm, 
w = 1.905 cm and c n = c(l + nj), n = 1, . . . , 15 with c = 2027.3 m/s and 7 = 0.02786. The 
width, height and length of the small cuboids are w' = 5.0 mm, h' = 5.0 mm, and e' = 6.0 mm, 
respectively. 

rods with l n = 1/(1 + wy), n = 1, . . . ,N and c n = (G/p) 1 ^ 2 were p is the density, G the 
shear modulus and I a fixed arbitrary length. Notice that, in circular rods, the velocity does 
not depend on the radius (l9| . For system B we take l n = I and c n = c(l + 717), c being an 
arbitrary constant velocity. The parameter 7 is dimensionless. 

To construct system B we use the Navier expression for the torsional velocity of cuboid 

n: 

(2) 



pin 

where /„ = (h n w 3 + h n w)/12 is the moment of inertia with respect to the axis of the system 
and a n is given by 
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For beams of varying h n we have verified Eqs. (J2J) and (j3J) experimentally. Solving these 
equations we have obtained the values of h n such that c n = c(l + nj). 
Then 



An) 

h 



( y/G/p(l + nj)j/2l for system A (4a) 
I c(l + nj)j/2l for system B, (4b) 



A, = Afj n) 



and the differences A/. = /j n+1 ' ) — /j n ^ are equal to 

\jGj ' frfj /2l for system A (5a) 
cyj/2l for system B, (5b) 

which are independent of index n. 

We shall now discuss system A. When the arbitrary parameter 7 is set equal to zero, 
a locally periodic rod is formed. This locally periodic rod shows a band spectrum {^(J. 
When 7 ^ a completely different spectrum occurs. The new spectrum then resembles the 
Wannier-Stark ladder. 

Before presenting the calculation of the normal modes for this system, and then showing 
numerical and experimental results, let us make a qualitative analysis to see what type of 
spectrum could be expected from the independent rod model. At the lowest frequencies, 
the wavelength A is of the same order of magnitude as L « J2n=i anc ^ the whole rod is 
excited. But when A decreases and becomes of the order of l\ — 1/(1 + 7), the longest rod, 
say rod 1, is excited in a state equivalent to its lowest normal mode. The rest of the A" 
rods are out of resonance, so the amplitude decreases as we move farther away from rod 1. 
Therefore, the state is localized around the latter. In some sense this was to be expected 
since we are disturbing a periodic structure to obtain a disordered one-dimensional system, 
which always shows localized wave amplitudes. Increasing the exciting frequency by A x 
the rod with length Z 2 = 1/(1 + 27), that is, rod 2, will now be excited and the rest will 
be out of resonance. The amplitudes of the vibrations therefore decrease as their distance 
from rod 2 increases. The wave amplitude is again localized but now around rod 2; it has a 
similar shape as the wave amplitude that rod 1 had before, but it has been slightly deformed, 
squeezed, and translated from rod 1 to rod 2. The same argument applies when rod n of 
length l n = 1/(1 + rvy) is excited. 

What we have done is to produce a finite WSL, i.e., A^ localized states with constant 
difference in frequency given by Eq. (|5ajl . However, more ladders exist since normal modes 
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FIG. 2: Block diagram of the experimental setup. For system A both detector and exciter are 
EMATs, while for system B the exciter was a piezoelectric transducer. 



with two or more nodes can also be excited in each rod. For instance, taking j = 2 in 
Eq. (|4*aj) a second ladder is obtained. This ladder is different from the first one because the 
frequency difference is now twice the one of the lower WSL, as can be seen from Eq. (|5a|). 
The states are again localized and all have a similar shape although squeezed. A third ladder 
exists with A3 = 3Ai and so on for other values of j. A similar argument for system B 
shows also the existence of several WSL. The difference between the quantum-mechanical 
and elastic ladders is that in the latter the spacing between resonances is not the same for 
different ladders. 

We have calculated the eigenmode properties of the rods of Fig. ^ (a) and (b) with free- 
end boundary conditions using the transfer matrix method for torsional waves discussed in 



Rcf. 



20j . The normal mode frequencies and amplitudes were measured using the experi- 



mental set up described in Fig. El We use an electromagnetic acoustic transducer (EMAT) 
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FIG. 3: Normal-mode frequencies of System A as a function of the dimensionless parameter 7. 



which is very versatile and operates at low frequencies. This EMAT, which we have recently 



developed 2l|, can selectively excite or detect compressional, torsional or flexural vibra- 
tions. In the inset of Fig. |2]the EMAT has been installed to detect torsional vibrations; see 
Ref. Q- 

We show in Fig. El the spectrum of system A as a function of the dimensionless parameter 
7. As mentioned above, for 7 = a band spectrum appears, and as 7 grows the levels of each 
band separate to form the WSL. The normal-mode frequencies of the rod shown in Fig. ^ 
(a) are given in Fig. 0] for 7 = 0.091. We first note that the theoretical results coincide 
very well with the experimental ones. Furthermore, the qualitative treatment provides a 
rather good first approximation. One can see from this figure that the states form a set 
of Wannier-Stark ladders, as discussed before. The first band composed by the extended 
modes is not displayed in this graph. Notice that the frequencies in the extremes of each 
ladder do not have the same difference in frequency as those at the middle of the ladder. 
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FIG. 4: (color online) Normal-mode frequencies of System A yielding the elastic Wannier-Stark 
effect. For each value of j the left hand side column corresponds to the experimental values, 
the middle column to the numerical results obtained using the transfer matrix method and the 
right hand side shows the approximate results following from the independent rod model. In the 
calculation we used an effective value for r/R 2(j]. The uncertainty in the experimental values is 
less than 0.01%. In the insets the theoretical wave amplitudes are given for a state at the extreme 
of the ladder and another one in the center of it. 

This is due to a border effect in the wave amplitudes localized near the free ends. As shown 
in the insets of Fig. 01 the border amplitude lacks a portion of the wave amplitude that the 
states at the center of the ladder have. 

In Fig. El we show the comparison of theoretical and experimental wave amplitudes. These 
are localized around rod n. For example, in Fig. El (a) the sixth rod resonates and in Fig.|3](b) 
another state corresponding to the same ladder is localized around the tenth rod. Both have 
the same form although squeezed. Figures El (c) and El (d) show two states of the second 
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FIG. 5: (color online) Wave amplitudes of system A for 19 nodes (a), 23 nodes (b), 33 nodes (c), 
and 41 nodes (d). The continuous line corresponds to the transfer-matrix results and the dots to 
the measurements. The vertical lines along the rod axis indicate the position of the notches. 

ladder, with n = 6 and n = 11, respectively. Localization is again observed and, as was to 
be expected, the amplitudes now show two nodes in the resonating rods. Note the excellent 
agreement between theory and experiment after adjusting the height of the theoretical wave 
amplitude at only one point. 

As mentioned above, system B shows similar properties. We present in Fig. El as an 
example, two wave amplitudes for the first WSL of system B. In contrast with system A, 
these wave amplitudes have the same shape but translated and are not squeezed. Notice 
that the one-dimensional transfer matrix method fit the experimental wave amplitudes in 
spite of the fact that w and h n are not much smaller than /. 
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FIG. 6: (color online) Two wave amplitudes of system B, (a) one localized on the third rod with 
frequency / = 44.256 kHz and (b) other localized on the tenth rod with frequency / = 51.258 kHz. 
The double vertical lines along the beam axis indicate the position of the notches. 

To conclude, we have constructed an elastic analogue of Wannier-Stark ladders. In 
contrast with the optical analogue of Refs. we have observed the WSL directly. 

Furthermore, we measured for the first time the wave amplitudes, including phases, which 
show localization. We also observed higher Wannier-Stark ladders. The elastic Wannier- 
Stark ladders have potential applications in the design of systems with localized vibrations. 

This work was supported by DGAPA-UNAM under projects IN104400 and IN104903-3 
and by CONACyT Mexico, project 41024-F. 
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